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analysis/synthesis  sinusoidal  representations  and  transformations  were  developed  as  two 
interactive  tools  with  Graphical  User  Interface  (GUI)  using  MATLAB.  In  addition,  an 
interactive  tool  for  signal  frequency  component  editing  based  on  the  sinusoidal  model  is 
also  presented  in  this  thesis. 
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I.  INTRODUCTION 

A.         OVERVIEW 

Sinusoidal  representation  is  a  useful  model  for  speech  and  sound  analysis/synthesis. 
It  has  been  shown  that  the  synthetic  waveform  preserves  the  general  waveform  shape  and  is 
perceptually  indistinguishable  from  the  original  sound  [Refs.  1,2,3]-  However,  in  a  number 
of  applications  it  is  required  to  transform  a  sound  signal  to  a  different  waveform  which  is 
more  useful  than  the  original.  For  example,  in  time-scale  modification  of  speech,  the  rate  of 
articulation  may  be  slowed  down  in  order  to  make  degraded  speech  more  comprehensible. 
Alternatively,  the  sound  can  be  speeded  up  so  we  can  quickly  scan  a  passage  or  compress  it 
into  a  fixed  time  interval.  In  other  applications,  the  sound  is  compressed  or  expanded  in 
time  or  frequency.  For  instance,  in  music  synthesis  it  is  useful  to  change  the  length  or  pitch 
of  a  tone  without  changing  its  tonal  quality  or  timbre.  In  all  of  these  cases,  it  is  desired  to 
perform  sound  modification.  This  thesis  implements  a  fixed  time-scale  modification  and 
frequency  scaling  based  on  the  sinusoidal  model  [Ref.  4].  The  above  sound 
analysis/synthesis  sinusoidal  representations  and  transformations  have  been  developed  as 
two  interactive  tools  with  Graphical  User  Interface  (GUI)  using  MATLAB. 

Since  some  frequency  components  in  a  signal  are  redundant  (they  may  either 
correspond  to  the  noise,  or  carry  no  information),  we  may  not  want  to  include  them  when 
we  resynthesize  the  sound  signal.  In  other  cases,  only  a  portion  of  the  original  signal  needs 
to  be  regenerated,  or  a  small  part  of  the  signal  is  required  to  be  repeated  at  a  specific  time 
instant.  In  these  applications,  we  need  to  have  the  ability  to  edit  the  frequency  components 
of  the  synthetic  signal.  Thus,  an  interactive  tool  for  signal  editing  based  on  the  sinusoidal 
model  is  also  presented  in  this  thesis. 


B.         THESIS  OUTLINE 

The  remainder  of  this  thesis  is  organized  as  follows.  Chapter  II  addresses  the  sine- 
wave  speech  model  in  two  parts.  First,  the  speech  analysis/synthesis  model  based  on  a 
sinusoidal  representation  is  presented  [Ref.  1].  Following  this,  algorithms  for  fixed  time- 
scale  modification  and  frequency  scaling  based  on  the  sinusoidal  representation  are 
introduced  [Ref.  4].  Chapter  III  describes  the  implementation  of  the  three  interactive  tools 
for  sound  analysis/synthesis  with  GUI.  The  methods  implemented  include  signal  editing, 
fixed  time-scale  modification,  and  frequency  scaling.  Results  of  their  use  are  also  shown 
here.  Finally,  Chapter  IV  gives  conclusions  and  recommendations  for  future  work. 


II.  ANALYSIS  /  SYNTHESIS  BASED  ON  A  SINUSOIDAL  REPRESENTATION 

A.         SINE-WAVE  ANALYSIS  /  SYNTHESIS  MODEL 

1.  The  Sinusoidal  Representation 

A  speech  modeling  technique  developed  by  McAulay  and  Quatieri  [Ref.  1],  is  based 
upon  a  sinusoidal  representation  of  the  original  waveform.  In  the  general  speech  production 
model,  the  output  waveform  is  assumed  to  be  the  result  of  passing  the  glottal  excitation 
e(t)  through  a  linear  time-varying  filter  h(t,x) which  models  the  vocal  tract.  The  model 
can  be  written  as 


s(t)=  J  e(x)h(t,t-x)dx  . 


(1) 


and  is  depicted  in  Figure  1. 


it) 


Figure  1.  Model  of  Speech  Production 


In  the  sinusoidal  model  the  excitation  is  written  as  a  sum  of  sinusoids  with  time-varying 
amplitudes  and  phases,  namely 


where  Qf(t)  is  given  by 


with 


Lt\ 


eit)  =  TaMC0S[QM]  ' 


n,M  =  xF,M+4>, 


(2) 


(3) 


(4) 


In  this  model,  tf  is  the  onset  time  of  the  t1  sine  wave  and  L(t)  is  the  number  of  sine-wave 

components  at  time    t .   For  the  f  sine-wave  component,    ae(t)  is  the  time- varying 

amplitude,  ©,,(/)  is  the  frequency  and  Qe(t)  the  phase  corresponding  to  the  f1  sine  wave. 

The  quantity  4^(0   is  the  time- varying  contribution  to  the  phase  while  <j>,  is  a  fixed  phase 

offset  needed  since  the  sine  wave  components  for  different  indices  I  are  generally  not 
aligned. 

If  the  vocal  tract  transfer  function  is  written  as 

H(a,t)  =  M(cM)exp[yO(cM)]  ,  (5) 

then  the  output  speech  s(t)  can  be  written  as 

L{t) 

where 

At(t)  =  ae(t)Mt(t)  ,  .       .  (7) 

and 

Qe(t)  =  Q,(t)  +  <D,(f)  ,  (8) 

are  the  amplitude  and  phase  of  the  &  sine  wave  component  corresponding  to  the  frequency 
co,(f). 

The  sine  wave  model  has  been  found  to  be  useful  for  modeling  other  types  of 
sounds  besides  speech.  Other  specific  applications  for  this  method  have  been  in  music 
synthesis  [Ref.  1,  2]  and  in  underwater  acoustics  [Ref.  3].  For  these  applications,  the 
separation  of  the  sound  into  excitation  and  system  components  as  shown  in  Figure  1  may  or 
may  not  be  appropriate.  Still,  most  of  the  basic  ingredients  of  the  model  remain  and  can  be 
applied  in  these  applications. 

2.  Analysis 

The  purpose  of  the  analysis  step  is  to  estimate  the  composite  amplitudes, 
frequencies,  and  phases  of  the  sine  wave  model.  This  can  be  done  from  the  high-resolution 


short-time  Fourier  transform  (STFT).  The  original  analysis  method  proposed  by  McAuley 
and  Quatieri  [Ref.  1]  uses  a  purely  sine-wave-based  model  (i.e.,  the  excitation  and  system 
contributions  of  each  sine-wave  component  are  not  explicitly  represented).  In  their 
following  work,  a  new  analysis  procedure  is  developed  which  separates  the  vocal  cord 
excitation  and  vocal  tract  system  contributions  as  described  above.  Since  we  are  interested 
in  more  general  types  of  sounds  rather  than  speech,  the  original  analysis  method  along  with 
the  modified  amplitude  and  phase  representations  is  used.  Thus,  we  account  only  for  the 
model  of  the  vocal  cord  excitation  contribution  and  ignore  the  vocal  tract  system 
contribution.  With  this  simplification,  Eq.  (7)  and  Eq.  (8)  become 

At(t)  =  ae(t)t  (9) 

and 

ee(t)  =  Qe(t)  .  (io) 

The  analysis  proceeds  as  follows.  First,  the  data  is  sectioned  into  frames  of  equal 
length  for- spectral  analysis  and  a  Hamming  window  is  applied  before  taking  the  Fourier 
transform.  Frames  are  formed  at  an  interval  of  less  than  the  frame  length  allowing  for 
overlap  of  data.  For  a  speech  signal,  a  frame  length  of  20  ~  30  milliseconds  and  overlap 
interval  of  10  ~  15  milliseconds  are  recommended  [Ref.  1].  If  the  Fourier  transform  of  the 

windowed  speech  segment  is  written  as  S((£>,kR),  then  the  frequencies  of  e{t)  in  Eq.  (2)  at 

time  kR  (i.e.,  the  A*  analysis  frame),  are  chosen  to  correspond  to  the  L(kR)  largest  peaks  in 

the  magnitude  of  the  short-time  Fourier  transform,  |jS(g>,£R)|.  The  locations  of  the  largest 
peaks  are  estimated  by  looking  for  a  change  of  slope  from  positive  to  negative  of  the 
Fourier  transform  magnitude. 

If  we  denote  the  frequency  estimate  of  the  Ith  sinusoidal  component  at  the  A* 
analysis  frame  by  co*  =co^(kR),  then  the  amplitudes  and  phases  of  the  sine-wave 
component  are  given  by  the  samples  of  S(cd,kR)  at  the  specific  frequency  positions.  In 
other  words,  the  amplitudes  and  phases  are  written  as 


akt=\s{Skt,kR) 


(11) 


and 


Qke=ar^s(3[,kR)]  ,  (12) 

where  "arg"  denotes  the  principal  phase  value.  A  block  diagram  of  the  analysis  scheme  is 
given  in  Figure  2. 
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Figure  2.  Block  Diagram  of  Sinusoidal  Analysis 

The  number  of  peaks  are  not  constant  from  frame  to  frame  in  general,  and  there  will 
be  spurious  peaks  due  to  the  effects  of  window  sidelobe  interaction.  In  addition,  the 
locations  of  the  peaks  will  change  as  the  pitch  changes;  and  rapid  changes  in  both  the 
location  and  number  of  peaks  often  occur  in  certain  regions  of  the  sound  signal.  In  order  to 
account  for  such  movements  in  the  spectral  peaks,  the  concept  of  "birth"  and  "death"  of 
sinusoidal  components  is  introduced  here.  Suppose  that  the  peaks  up  to  frame  k  have  been 
matched  and  a  new  parameter  set  for  frame  k  +  1  is  generated.  We  now  attempt  to  match 

frequency  co*  in  frame  k  to  the  frequencies  in  frame  k  +  1.  If  all  frequencies  co*+1  in  frame 

k  +  1  lie  outside  a  "matching  interval"  of  co* ,  then  the  frequency  track  associated  with  co* 
is  declared  "dead"  on  entering  frame  k  +  1 .  When  all  frequencies  of  frame  k  have  been 
tested  and  assigned  to  continuing  tracks  or  to  dying  tracks,  there  may  remain  frequencies  in 

frame  k  +  1  for  which  no  matches  have  been  made.  It  is  assumed  that  such  frequencies  co  *„+1 

were  "born"  in  frame  k  and  a  new  frequency  co  *  is  created  in  frame  k  with  zero  magnitude. 


This  procedure  is  done  for  all  unmatched  frequencies.  Further  details  of  this  "birth"  and 
"death"  matching  procedure  can  be  found  in  [Ref.  1]. 

The  result  of  applying  this  method  to  a  segment  of  a  sound  signal  is  shown  in  Figure 
3.  Each  horizontal  line  represents  a  particular  frequency  component  which  is  present  for 
some  number  of  frames.  These  lines  are  called  "frequency  tracks."  The  frequency  tracks 
demonstrate  the  ability  of  the  method  to  adapt  quickly  through  the  transitory  regions  such  as 
voiced/unvoiced  transitions  in  speech.  Typically  there  are  many  very  short  frequency  tracks. 
Some  of  these  may  not  contribute  significantly  to  the  general  structure  of  the  waveform  but 
merely  serve  to  match  small  details.  As  will  be  seen  later,  the  editing  tools  developed  in  this 
thesis  allow  one  to  eliminate  many  of  these  shorter  frequency  tracks  and  thus  simplify  the 
sinusoidal  model  for  the  signal. 


4000 
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Figure  3.  Typical  Frequency  Tracks  for  a  Sound  Signal 


3.  Synthesis 

Sound  signal  synthesis  from  the  sine-wave  parameters  begins  with  matching  the 
amplitude  and  phase  samples  in  Eq.  (11)  and  Eq.  (12)  of  each  sine-wave  computed  at 
consecutive  frame  boundaries.  This  is  followed  by  interpolation  of  the  resulting  pairs  of 
amplitude  and  phase  samples  of  the  signal  over  each  frame.  The  interpolation  of  parameters 
is  based  on  the  assumption  that  the  signal  is  "slowly  varying"  across  each  frame  and  that  the 
frequencies  of  the  sine  waves  form  smooth  frequency  tracks  co  t  (r) .  This  constraint  allows 
us  to  interpolate  samples  over  a  frame  duration.  If  linear  interpolation  is  used  for  the 
amplitude,  the  amplitude  estimate  ae{i)  over  the  A*  frame  is  given  by 

ae(t)  =  a$+(a™-ake)t/T  ,  (13) 

where  a)  and  <S*+1   are  a  successive  pair  of  excitation  amplitude  estimates  for  the  & 

frequency  track,  T  is  the  frame  duration  and  t  e  [0,  T]  is  the  time  into  the  k^  frame. 

This  simple  linear  interpolating  procedure  cannot  be  used  for  estimating  the  phase 

and  frequency  of  the  sinusoid  over  a  frame,  however.  This  is  because  the  phase  Q*  may 

contain  discontinuities  of  In  since  the  phase  of  S((£>,kR)  in  Eq.  (12)  is- measured  modulo 

2n.  Hence,  phase  unwrapping  must  be  performed  for  interpolation  of  the  excitation  phase  to 
ensure  that  the  frequency  tracks  are  sufficiently  "smooth"  across  the  frame  boundaries.  A 
cubic  polynomial  for  solving  this  problem  is  first  proposed  in  [Ref.  1]  for  sine- wave-based 
synthesis.  For  the  duration  at  a  single  frame  the  estimate  is  defined  as 

Cle(t)  =  a  +  bt  +  ct2  +  dt3  ,  (14) 

with  t  =  0  corresponding  to  frame  k  and  t  =  T  corresponding  to  frame  k  + 1 .  The 
instantaneous  frequency  is  then  the  derivative  of  the  phase,  namely 

© ,(/)  =  —  &({t)  =  b  +  2ct  +  3dt2  .  (15) 

at 

In  order  to  provide  a  good  synthetic  waveform,  it  is  necessary  that  the  cubic  phase  function 
and  its  derivative  equal  the  excitation  phase  and  frequencies  measured  at  the  frame 


boundaries.  By  using  the  algorithms  in  [Ref.  1],  the  resulting  phase  function  not  only 
matches  the  phase  at  the  frame  boundaries,  but  also  resolves  the  In  phase  discontinuities. 
Details  of  phase  unwrapping  and  cubic  interpolation  can  be  found  in  [Ref.  1]. 

It  was  noted  earlier  that  the  phase  estimate  over  the  \t  frame  can  be  written  in  terms 
of  a  time-varying  term  and  a  constant.  Specifically,  from  Eq.  (3)  and  Eq.  (4) 

Cle(t)=  jo^GJda  +  h 

=  j6e(a)da+  J  a>^(a)dcr  +  <j)£  , 

where  the  time  origin  (t  =  0)  is  taken  to  be  at  the  beginning  of  the  current  frame  and  te  is 

the  onset  time  of  the  £*  sine-wave.  Let  V  denote  the  phase  due  to  the  time-varying 
frequency  accumulated  up  to  frame  t,  that  is, 

Z*    =\t®Md<*  ■  (17) 

If  Vt(t)  denotes  the  phase  due  to  the  time-varying  frequency  accumulated  over  frame  kr, 
that  is, 

^M=£©*M<fe>  as) 

then  the  excitation  phase  can  be  written  as 

6,W  =  ^W+Z*+$<  •  (i9) 

The  resulting  excitation  phase  function  Qe(t)  consists  of  a  constant  component  and 
a  time-varying  portion.  The  constant  component  consists  of  two  parts:  the  phase  offset 

estimate  tyf,  and  the  accumulated  phase  component  ^  ,  which  can  be  obtained 
recursively  as 

The  interpolated  amplitudes  and  phases  are  used  to  generate  sinusoids  which  are  then 
summed  to  generate  the  output  sound  signal.  The  final  synthetic  waveform  is  written  as 


L(t) 


^(0=E^(Ocos[n,(0] 


where 


A  block  diagram  of  the  synthesis  structure  is  given  in  Figure  4. 
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Figure  4.  Block  Diagram  of  Sinusoidal  Synthesis 

B.         TIME-SCALE  AND  FREQUENCY  TRANSFORMATION 

1.  Fixed  Rate  Change 

The  goal  of  time-scale  modification  is  to  maintain  the  perceptual  quality  of  the 
original  sound  while  changing  the  apparent  rate  of  sound  production.  In  speech,  the 
technique  is  used  to  synthesize  speech  corresponding  to  a  person  speaking  more  rapidly  or 
slowly  without  changing  the  quality  of  the  person's  voice.  The  scheme  illustrated  here  is 
based  upon  the  algorithm  developed  by  Quatieri  and  McAulay  with  slight  simplification 
[Ref.  4].  Although  the  authors  proposed  both  fixed  rate  change  and  time-varying  rate 
change,  only  the  fixed  rate  change  is  performed  here. 
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For  a  fixed  time-scale  transformation,  the  time  t0  corresponding  to  the  original 
sound  production  rate  is  mapped  to  the  transformed  time  t'Q  through  the  mapping  t'Q  =  pr0 . 
The  case  p  >  1  corresponds  to  time-scale  expansion,  while  the  case  p  <  1  corresponds  to 
time-  scale  compression.  The  case  of  time-scale  expansion  is  depicted  in  Figure  5. 

U 


'0 

*~t 


0 

P>1 


Jl — *f 


0  4=P'o 

Figure  5.  Time  Warping  with  Fixed  Rate  Change  p  >  1 

In  the  sine-wave  model  discussed  here,  the  parameters  which  are  scaled  are  the 
model  amplitudes,  frequencies,  and  phases.  The  model  parameters  are  modified  so  that 
frequency  tracks  (0  ((t)  are  stretched  or  compressed  in  time  while  the  value  of  (0((t), 
which  corresponds  to  pitch,  is  maintained.  The  mathematical  model  for  the  fixed  time-scale 
modified  sound  s'(t') ,  is  then  given  by 

s'{t')=Yja[{t')cos[Q.\{t')]  ,  (23) 

where 

a'l(t,)  =  al(p-lt')  ,  (24) 

and 

Cl'e(t')=  jt(ae(p~XT)dx  +  ^i  .  (25) 

Letting  a  =  p_Ix  ,  then  Q'e (/')  can  be  written  as 

Q;(r)=  f      (De(<j)da/p~]  +<|>, 

(26) 

=  r«(p-V)/p-+(z;)  +♦, . 
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Since  these  model  parameters  are  derived  on  a  frame-by-frame  basis,  we  can  think  of  the 
inverted  time  p~V  as  the  time  into  the  A*  frame  within  the  original  time  scale.  Therefore, 
the  fixed  time-scale  synthetic  waveform  can  be  obtained  as 


L(t') 


where 


and 


£=1 


a[{t')  =  a\{9-xt')T 


(27) 


(28) 


Q.'Xt')  =  v[{9-lt'). 
and  I  ^     J   computed  recursively  as 


P-+IZ 


+* 


a   •> 


(29) 


(30) 


The  notation  (  )T  in  Eq.  (28)  denotes  modulo  T,  which  is  the  original  frame  duration.  Figure 
6  illustrates  an  example  in  which  a  segment  of  male  speech  is  expanded  by  a  factor  of  2. 

2.  Frequency  Scaling 

The  sound  can  be  changed  in  pitch  by  performing  frequency  scaling.  This  is 
accomplished  by  taking  the  synthesized  phase  to  be 

where  P  is  the  scaling  factor  for  each  frequency  track  co  t  (f ) .  The  operation  performed  here 
is  equivalent  to  shifting  the  frequency  tracks  to  new  locations.  The  resulting  modified 
waveform  over  the  ^  frame  is  given  by 


(31) 


L{t) 


S'M  =  Z^Mcos[6;(/)] 


(32) 
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where 


6;to  =  pr,M+(X;*  )  +$,  ,  (33) 

with  I  ^        computed  recursively  as 


This  waveform  modification  corresponds  to  an  expansion  or  compression  of  frequency  and 
a  change  in  pitch.  Figure  7  illustrates  an  example  in  which  the  pitch  of  a  male  speech  is 
scaled  by  a  factor  of  2. 


13 


M 

i 

03 

| 

1 

i 

| 

0.2 

J  .    .. 

iT 

01 

In 

IS 

I  1 

Hi 

ii' 

II 

0 

i 

V'f  : 

M' 

jfcj! 

fro 
f 

'^1|i: 

0  1  - 

i 

1 

ftf14ll 

i 

0.2- 

i  j 

4 

:.  i 

03- 

i 

i 

1 

062  0.64  0.66 


(a) 


1  .    , 

1 

i 

If  til 

; 

1 

1 1  '■■■ 

11 

1 

1 

i,  1 

1 

I'l 

1 

I; 

[i 

1 

(fl 
if 

| 

1 

fl/i 

i 
i 

ill 

F  ' 

fl 

1' 

1 

i 

1 ! 
i  i'i 

[1 
tJ 

''-,  i 

I"" 

In 

ft 

!  ( 

1 

1 

■ 

1  1 

1 

i    ! 

1 

| 

t 

i 

(b) 

Figure  6.  Time-scale  Expansion  of  Speech  (a)  Original  (b)  Expansion  (  p=2) 


14 


I  :"■  ' 


■l: 


i 


% 


■!*■■'! 


U. 


i 


1.01  1.02  1.03 


1.04 
Time  (sec) 


1.05  1.0B  1.07  1.08 


(a) 


i 

i 

i 

i         i 

i 

i  . 

| 

li 

| 

i 

If  if 

ii 

i 

;  1 
fllj 

■Hi 

I 

11 

Ill 

Ijl 

;  : 

ii   | 

lii 

II 

M   1 

|il') 

i 

r  1    '' 

1 

1 

! 

11  ' 

All: 
in") 

1 
iff 

1   It 
1!  Ill 

■f  1 

P 

I! 

I 

■# 

r  i 

i> 

.[Ifl 

i 

1 

|| " 

'     ji 

j  I 

1 

* 

1 

! 

.   . 

i 

i         i 

i 

I 

- 

1.02  1.03  1.04  1.05  1.06  1.07  1.1 

Time  (sec) 


(b) 

Figure  7.  Frequency  Scaling  of  Speech  (a)  Original  (b)  Pitch-scaled  (  p=2) 


15 


16 


III.    IMPLEMENTATION  OF  THE  INTERACTIVE  TOOLS 

The  definition  of  a  user  interface  is  moving  from  a  command-line  oriented  interface 
to  one  that  includes  graphic  features.  Over  the  years,  graphical  user  interfaces  (GUI)  have 
grown  in  popularity.  GUI  use  push  buttons,  editable  boxes,  and  other  graphical  controls 
which  can  be  activated  with  a  mouse  to  select  various  options  and  execute  commands 
[Ref.  5].  The  purpose  here  was  to  develop  interactive  user  interface  tools  which  can  be  used 
to  perform  the  sound  analysis/synthesis,  frequency  track  editing,  and  sound  transformations 
based  on  the  sinusoidal  representation.  The  use  of  these  GUI  relieves  the  user  of  the  need  to 
memorize  a  large  number  of  textual  commands,  and  allows  liim/her  to  see  the  results  almost 
immediately.  The  interactive  tools  described  here  were  developed  on  Unix  workstations, 
and  require  MATLAB  version  4.2c  as  well  as  its  Signal  Processing  Toolbox.  Although 
MATLAB  provides  the  necessary  support  for  the  GUI  on  both  Unix  and  IBM  PC- 
compatible  platforms,  some  modifications  will  need  to  be  made  if  the  user  wants  to  use 
these  tools  on  IBM  compatible  PCs. 

A.         THE  SOUND  ANALYSIS  /  SYNTHESIS  INTERACTIVE  TOOL 

An  interactive  sound  analysis/synthesis  tool  based  on  sinusoidal  representation 
model  is  described  in  this  section.  This  tool  allows  the  user  to  analyze  an  existing  sound 
waveform,  extract  the  parameters  that  represent  a  quasi-stationary  portion  of  that  waveform, 
and  then  use  those  parameters  to  reconstruct  an  approximation  that  is  "very  close"  to  the 
original  signal.  In  other  words,  the  algorithms  behind  this  tool  contain  two  parts,  namely, 
analysis  and  synthesis.  When  this  tool  is  invoked,  the  user  must  indicate  a  sound  signal  as 
an  input  argument  in  the  associated  .m  function;  the  signal  will  be  drawn  in  the  top  portion 
of  the  window  as  shown  in  Figure  8  (a)  and  labeled  "Original  Signal."  After  loading  the 
signal  into  the  workspace,  the  user  needs  to  provide  some  important  values  which  are  used 
in  the  signal  analysis  and  synthesis  algorithms. 
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The  first  value  is  the  sampling  frequency  for  the  analysis  and  synthesis  procedures 
which  should  be  same  as  the  value  used  in  digitizing  the  original  sound  signal.  For  all  cases 
discussed  in  this  thesis,  the  sampling  frequency  8,000  Hz  is  used.  This  is  close  to  the  actual 
value  of  8,192  Hz  used  in  the  SUN  Unix  workstations. 

The  next  values  to  be  entered  are  the  windowed  frame  length  and  overlap  width.  A 
windowed  frame  greater  than  20  milliseconds  is  sufficient  for  generating  a  good  quality 
synthetic  waveform  according  to  [Ref.  1];  this  corresponds  to  160  points  if  the  sampling 
frequency  is  8,000  Hz.  A  50%  overlap  of  the  frame  is  recommended  for  this  sinusoidal 
representation  model,  which  would  result  in  a  frame  overlap  width  of  80  points.  The  default 
values  for  the  windowed  frame  length  and  overlap  width  in  this  tool  are  200  and  1 00  points, 
respectively.  These  two  user-input  values  are  shown  in  the  second  and  third  editable  boxes 
in  Figure  8  (a)  and  (b). 

The  fourth  parameter  is  the  threshold  level  (in  dB),  which  allows  the  user  to  limit 
the  maximum  number  of  peaks  detected  over  a  frame.  The  typical  range  for  this  threshold 
value  is  from  60  to  90  dB.  A  default  value  of  80  dB  has  been  used  throughout  the 
experiments.  In  general,  the  performance  will  not  be  affected  much  by  the  choice  of  this 
threshold  level  unless  too  few  peaks  are  allowed. 

A  concept  of  "birth"  and  "death"  of  sinusoidal  components  was  described  earlier  in 
Chapter  I  (B)  and  is  used  to  account  for  the  rapid  change  on  both  the  number  and  location 
of  spectral  peaks.  The  fifth  input  value  in  this  tool  is  the  frequency  interval  used  while  the 
frame-to-frame  peak  matching  procedure  is  performed.  It  indicates  the  number  of  frequency 
bins  that  the  frequencies  on  two  successive  frames  can  deviate  and  still  be  considered  to  be 
"matched."  A  value  of  1 0  has  been  set  as  a  default. 

The  last  input  value  is  the  number  of  points  used  in  the  computation  for  discrete 
Fourier  transform  (DFT)  of  each  frame.  Typically,  512  to  1024  points  should  be  enough  for 
generating  the  synthesis  signal  if  the  frame  length  does  not  exceed  500  points.  In  this  thesis, 
all  experiments  were  done  using  1024  points. 
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Having  entered  these  values  the  user  now  is  ready  to  do  the  sound  signal 
analysis/synthesis.  The  "Synthesize"  push  button  activates  both  the  analysis  and  the 
synthesis  functions.  In  other  words,  when  the  user  presses  the  "Synthesize"  button,  the  tool 
extracts  the  waveform  parameters  first,  which  are  required  for  the  model,  and  then  passes 
those  parameters  to  the  synthesis  system  so  that  the  synthetic  waveform  will  be  generated. 
The  result  after  entering  the  parameters  and  pressing  the  "Synthesize"  button  is  as  shown  on 
the  bottom  portion  in  Figure  8  (b). 

Of  the  above  six  user-input  values,  only  the  first  three  values  (i.e.,  sampling 
frequency,  windowed  frame  length,  and  overlap  width)  are  essential  and  case-dependent  for 
the  sound  signal  analysis/synthesis.  It  is  usually  not  necessary  to  change  the  other  three 
values  (i.e.,  threshold  level,  frequency  matching  interval,  and  DFT  points).  Additionally, 
push  buttons  are  available  for  the  users  to  "play"  (i.e.,  listen  to)  both  the  original  and 
synthesis  sound  signal  using  the  platform's  audio  output.  For  new  users  and  users  that  are 
not  familiar  with  this  interactive  tool,  an  on-line  help  function  is  available  by  pressing  the 
"Help"  push  button. 

B.         THE  FREQUENCY  TRACK  EDITING  INTERACTIVE  TOOL 

The  frequency  track  editing  tool  allows  the  user  to  "edit"  frequency  components  of  a 
signal  by  inputting  some  appropriate  parameters  and  even  by  using  a  pointing  device  such 
as  a  mouse.  In  a  number  of  applications,  it  is  desired  to  make  the  synthetic  signal  fit  in  a 
specific  time  interval,  or  to  eliminate  short  frequency  tracks  to  simplify  the  model. 
Frequently,  some  of  these  short  tracks  only  serve  to  match  the  "detail"  in  the  original 
waveform,  and  removing  them  will  not  change  the  general  characteristics  of  the  sound. 

Since  the  results  of  the  sinusoidal  representation  model  are  very  robust,  it  is  not 
necessary  to  reconstruct  the  signal  with  all  frequency  components  which  are  extracted  from 
the  original  waveform.  This  interactive  tool  offers  five  frequency  track  editing  functions  for 
users  who  wish  to  generate  different  types  of  synthetic  signals. 
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The  first  option  provided  by  this  tool  allows  the  user  to  eliminate  all  frequency 
tracks  of  less  than  some  specified  length.  After  eliminating  these  tracks,  the  associated 
signal  is  resynthesized  according  to  the  new  frequency  tracks.  Figure  9  illustrates  an 
example  in  which  all  frequency  tracks  less  than  20  frames  in  length  are  removed.  Notice 
that  there  is  no  large  difference  between  the  "original"  synthesis  waveform  in  (a)  and  the 
"new"  synthesis  waveform  in  (b),  although  the  underlying  model  has  been  considerably 
simplified. 

In  many  signal  processing  applications,  users  are  likely  to  design  different  kind  of 
filters,  such  as  low  pass,  high  pass,  or  band  pass  filters,  in  order  to  eliminate  the  unwanted 
frequency  components  and  preserve  the  specific  range  of  frequencies  which  carries  the 
information  they  need.  The  second  option  offered  by  this  tool  allows  the  user  to  implement 
those  filters  very  easily,  just  by  indicating  the  specific  range  of  frequencies  to  be  removed. 
The  example  in  Figure  1 0  illustrates  the  result  of  eliminating  frequencies  in  the  range  from 
3,000  Hz  to  4,000  Hz,  which  corresponds  to  passing  the  signal  through  a  low-pass  filter, 
and  the  associated  "new"  synthetic  waveform.  In  this  case  there  are  not  many  long  tracks  of 
high  frequency  components  in  the  range  of  frequencies  removed  so  there  is  only  a  slightly 
noticeable  effect  on  the  waveform.  The  elimination  of  the  higher  frequencies  are  most 
apparent  when  listening  to,  or  "playing"  the  sound. 

Another  example  is  shown  in  Figure  1 1  (b)  where  all  frequency  tracks  less  than  20 
frames  in  length  are  removed,  followed  by  eliminating  the  frequency  range  from  3,000  Hz 
to  4,000  Hz.  Figure  11  (a)  again  shows  the  original  frequency  track  plot  and  associated 
synthetic  waveform. 

In  some  cases,  it  is  desired  to  regenerate  the  signal  using  only  part  of  the  original 
signal.  Thus,  we  may  need  to  be  able  to  "cut"  a  small  region  in  time  of  the  frequency  tracks, 
and  then  regenerate  the  signal  again.  This  tool  allows  the  user  to  indicate  the  frame  range  to 
be  cropped.  An  example  is  shown  in  Figure  12  where  frames  30  to  40  have  been  cut, 
therefore  eliminating  some  of  the  "silent  region"  between  the  two  major  portions  of  the 
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sound.  The  resynthesized  waveform  is  also  shown  in  the  top  portion  of  the  window  in 
Figure  12  (b). 

In  other  situations,  the  user  may  desire  to  move  or  repeat  a  small  portion  of  sound 
signal  at  a  specific  time  instant.  This  can  also  be  done  by  using  the  frequency  track  editing 
tool.  Users  are  required  to  input  the  range  of  frames  to  be  repeated,  and  the  position  in  time 
where  the  frames  are  to  be  inserted.  Figure  1 3  illustrates  an  example  in  which  frames  from 
30  to  40  are  copied  and  re-inserted  at  frame  30.  In  this  case,  the  result  is  an  increase  the 
"silent  region"  between  the  two  major  portions  of  the  sound.  The  resynthesized  signal  is 
shown  in  the  top  portion  of  the  window  in  Figure  13  (b). 

Sometimes,  it  is  desired  to  remove  some  specific  longer  frequency  tracks  after  most 
of  the  shorter  tracks  have  been  removed.  Many  times  this  is  not  possible  using  the  methods 
that  have  been  previously  described.  A  handy  mouse  frequency  track  editing  function  was 
developed  for  this  purpose.  The  user  activates  this  editing  function  by  pressing  the 
"START"  push  button  at  the  right  bottom  corner  of  the  window.  The  cursor  changes  from 
an  arrow  to  cross-hairs  indicating  that  the  editing  function  is  active.  The  user  then  places  the 
cross-hair  cursor  on  a  frequency  track  to  be  eliminated  and  "selects"  that  track  with  the  left 
mouse  button.  Every  selected  frequency  track  changes  from  its  normal  yellow  solid  color  to 
a  red,  dotted  line.  The  unaffected  frequency  components  will  be  saved  for  the  use  of 
generating  the  new  sound  signal. 

Users  are  allowed  to  select  as  many  frequency  tracks  as  they  wish  to  eliminate.  The 
"new  modified"  synthetic  signal  is  not  generated  until  the  user  has  finished  the  frequency 
track  selecting  step.  When  the  user  has  finished  selecting  frequency  components  he/she 
presses  the  right  mouse  button  in  the  region  and  then  presses  the  "OK"  button  to  synthesize 
the  waveform.  An  example  of  the  use  of  the  tool  and  this  function  is  shown  in  Figure  14.  In 
Figure  14  (a),  all  frequency  tracks  less  than  20  frames  in  length  were  initially  removed. 
Figure  14  (b)  shows  the  result  where  three  specific  additional  frequency  tracks,  one  at 
approximately  3,700  Hz,  one  at  2,300  Hz  and  one  near  0  Hz,  have  been  eliminated.  The 
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Figure  13.  Frequency  Editing  Tool  Sample  View  (a)  Original  (b)  After  Frequency 
Tracks  Frames  From  30  to  40  Have  Been  Repeated  at  Frame  30 
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tracks  eliminated  are  shown  as  dotted  lines  and  the  waveform  shown  in  Figure  14  (b)  is  the 
result  of  synthesis  of  the  removal  of  these  tracks. 

C.         TIME  AND  FREQUENCY  SCALING  INTERACTIVE  TOOL 

The  time  and  frequency  scaling  tool  provides  a  means  of  generating  the  expansion 
or  compression  of  a  sound  signal  in  the  time  domain,  as  well  as  a  change  in  spectral 
envelope  and  pitch  contour  according  to  the  methods  described  in  Chapter  II  (B).  There  are 
two  options  offered  by  this  tool.  The  first  is  time-scale  modification.  In  the  case  of  time- 
scale  modification,  the  new  sound  signal  is  expanded  or  compressed  depending  on  the  value 
which  the  user  inputs.  The  new  sound  signal  is  expanded  if  the  value  is  greater  than  1 ,  and 
is  compressed  if  the  value  is  less  than  1 .  The  modified  signal  is  automatically  generated 
right  after  the  user  inputs  the  time-scale  factor.  An  example  is  illustrated  in  the  top  portion 
of  the  window  as  shown  in  Figure  1 5  where  a  segment  of  male  speech  is  expanded  by  a 
factor  of  2.  In  this  case,  it  is  found  that  the  rate  of  articulation  has  been  slowed  down  while 
the  perceptual  quality  of  the  original  sound  is  maintained. 

If  the  user  wishes  to  perform  a  frequency  transformation,  then  the  second  option  of 
this  tool  can  be  invoked.  The  user  can  increase  the  pitch  of  the  synthesized  sound  signal  by 
entering  a  pitch-scaling  factor  which  is  greater  than  1  or  lower  the  pitch  by  inputting  a  value 
that  is  less  than  unity.  Both  the  time-  and  pitch-scaling  factors  have  been  limited  to  the 
values  in  range  from  0.1  to  2,  since  values  outside  of  this  range  generally  produce  poor 
results.  An  example  of  frequency  scaling  by  a  factor  of  2  for  male  speech  is  depicted  in  the 
bottom  portion  of  the  window  shown  in  Figure  1 6.  The  resulting  speech  sounds  like  a 
young  boy's  voice  since  the  pitches  of  children's  voices  are  higher  than  those  of  adults  in 
general. 
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IV.  CONCLUSIONS 

A.         DISCUSSION  OF  RESULTS 

In  this  thesis,  a  sinusoidal  representation  model  for  sound  signals  by  McAulay  and 
Quatieri  is  described  and  is  used  in  an  analysis/synthesis  technique  based  on  the  amplitudes, 
frequencies,  and  phases  of  excitation  contributions  of  the  sine  wave  components.  In  the 
analysis  steps,  the  data  is  first  sectioned  into  frames  and  the  discrete  Fourier  transform 
(DFT)  is  applied  over  each  frame.  The  peaks  in  the  resultant  spectrum  determine  the 
frequencies  of  sinusoids  to  be  used  in  the  model  and  which  are  "tracked"  through 
successive  frames.  The  amplitudes  and  phases  of  the  sinusoids  are  given  by  the  appropriate 
samples  of  the  DFT  corresponding  to  those  peak  frequencies. 

In  the  synthesis  step,  these  amplitude  and  phase  functions  are  applied  to  the  sine- 
wave  generator,  which  adds  all  sinusoidal  components  to  produce  the  synthetic  signal 
output.  We  have  found  that  this  model  reproduces  the  sound  very  accurately  and  confirms 
the  claim  by  the  authors  that  the  sound  is  "perceptually  indistinguishable"  from  the  original 
sound  [Ref.  1]. 

Functional  relationships  for  each  of  the  sine- wave  parameters  have  been  developed 
by  the  original  authors  that  allow  the  synthesis  system  to  perform  a  variety  of  sound  signal 
transformations,  such  as  time-scale  modification  and  frequency  scaling  implemented  in  this 
thesis  [Ref.  4].  These  were  also  found  to  be  effective. 

All  of  the  above  sound  analysis/synthesis  sinusoidal  representations  and 
transformations  were  developed  as  two  interactive  tools  with  GUI  using  MATLAB.  In 
addition,  an  interactive  tool  for  signal  frequency  track  editing  based  on  the  sinusoidal  model 
was  also  implemented  so  users  can  simplify  the  model  and  reduce  its  complexity.  Examples 
of  the  use  of  these  tools  are  presented  in  this  thesis. 
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B.         SUGGESTIONS  FOR  FUTURE  STUDY 

In  the  sine-wave-based  time  and  frequency  scaling  modification  system,  the 
modified  synthetic  waveforms  are  good  in  perceptual  quality  but  some  structure  of  the 
original  waveform  is  lost.  A  new  sine-wave-based  speech  modification  algorithm  called  the 
"Shape  Invariant"  technique  has  been  proposed  by  the  original  authors  which  is  able  to 
maintain  the  temporal  structure  of  the  original  waveform  [Ref.  6].  It  would  be  worthwhile 
to  implement  this  new  technique  and  incorporate  it  into  the  existing  GUI. 

Although  this  sinusoidal  representation  model  produces  very  accurate  results,  it  is 
demanding  in  terms  of  computation.  Another  worthwhile  endeavor  would  be  to  improve  the 
computational  performance  of  the  sine-wave-based  modification  system.  MATLAB  is  poor 
in  executing  "Loop  Iteration"  code  because  it  is  an  interpreted  language.  Unfortunately,  the 
implementation  of  the  modification  system  uses  several  "loops"  which  makes  the  execution 
time  quite  long.  Rewriting  the  code  to  improve  the  computational  efficiency  would  be  very 
desirable  so  that  the  user  can  see  the  results  more  quickly.  Perhaps  compiling  the  code  with 
the  MATLAB  to  C  (or  C++)  compiler  would  result  in  faster  execution. 
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APPENDIX 

This  appendix  includes  the  sound  signal  examples  and  main  programs  described  in 
this  thesis.  Several  music  and  speech  signals  were  used  in  the  analysis/synthesis,  frequency 
track  editing,  and  transformation  experiments,  only  two  of  them  are  presented  in  this  thesis, 
however.  They  are: 

•  The  speech  phrase  "baseball"  from  a  male  speaker  (file:  obase.mat),  and 

•  Two  notes  (file:  tbhi2.mat)  excerpted  from  a  four-note  trombone  passage 
(file:  tb_hi.au). 

The  program  is  written  entirely  in  MATLAB  and  makes  extensive  use  of  Graphical 
User  Interface  (GUI)  features  from  MATLAB.  In  addition,  the  MATLAB  Signal  Processing 
Toolbox  is  required  to  run  these  programs. 

The  programs  are  divided  into  three  parts,  the  sound  analysis/synthesis,  frequency 
track  editing,  and  sound  transformation  (including  time-scale  modification  and  frequency 
scaling).  They  correspond  to  Chapter  III  (A),  (B),  and  (C),  respectively. 

A.         SOUND  ANALYSIS/SYNTHESIS  FUNCTIONS 

[cand,  mag,  phas,  sig,  par,  synt]  =  gsinwave(signal); 

This  is  the  main  program  which  calls  other  associated  MATLAB  functions  in  order 
to  perform  the  sound  analysis/synthesis.  The  user  invokes  this  program.  This  function  calls 
guisynth.m  which  implements  all  GUI  actions  and  calls  the  two  other  functions  analy.m  and 
synth.m  which  perform  the  anlysis  and  synthesis,  respectively.  The  input  to  gsinwave.m  is 
the  original  sound  signal  and  the  outputs  are  the  synthesized  signal  (synt)  and  a  set  of 
variables  which  are  generated  from  the  analy.m  and  synth.m  functions.  These  other 
variables  are  described  below. 
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1.  Analysis  Function 


[cand,  mag,  phas,  sig,  par]  =  analy(signal,  sam_fs,  N,  n,  peakjhr,  mat_win,  fft_size); 


Inputs 


Outputs 


Default  Values 

signal 

original  sound  signal 

sam_fs 

sampling  frequency 

8000  Hz 

N 

windowed  frame  length 

200  points 

n 

overlap  width 

100  points 

peak_thr 

peak  picking  threshold 

80  dB 

mat_win 

frequency  matching  interval 

10 

fft_size 

DFT  points 

1024 

cand 

sinusoids  peak  candidate  matrix 

mag 

magnitude  matrix  of  DFT 

phas 

phase  matrix  of  DFT 

sig 

windowed  original  signal 

par 

parameters  used  by  other  functions 

The  output  variables  listed  above  can  be  described  in  more  detail  as  follows. 


•     cand 


This  matrix  contains  the  "matched"  (i.e.,  after  applying  the  "birth"  and  "death" 
process  and  "frame-to-frame"  peak  matching  procedure) "  peak  frequencies 
information  of  sinusoids  which  are  extracted  from  the  DFT  spectrum.  The 
number  of  rows  of  this  matrix  is  equal  to  one  half  of  the  DFT  points,  and  the 
number  of  columns  depends  on  both  the  windowed  frame  length  and  length  of 
the  sound  signal. 

A  cand  matrix  of  the  sound  tbhi2.mat  was  generated  by  using  the  default  input 
values  mentioned  earlier  in  the  analysis  function.  The  size  of  this  cand  matrix  is 
512  by  77.  A  small  matrix  corresponding  to  rows  36  to  60  and  columns  1  to  10 
was  excerpted  from  the  cand  matrix  as  an  example.  The  location  and  value  of 
each  element  indicate  the  frequency  and  the  "matched"  position  on  the 
following  frame,  respectively. 

For  instance,  the  element  (37,2)  of  the  original  cand  matrix  corresponds  to  the 
frequency  value  289  Hz.  The  value  "38"  in  this  position  of  the  matrix  indicates 
that  the  element  (38,3)  is  the  next  matched  position  ("38"  corresponds  to  the 
frequency  value  297  Hz).  Also,  since  there  is  no  value  "37"  in  column  1  of  this 
matrix,  this  frequency  track  is  considered  to  be  "born"  in  frame  1.  Let  us 
examine  another  element  (53,2)  of  this  matrix.  Its  value  "53"  indicates  that  the 


element  (53,3)  would  possible  be  a  nonzero  value  so  that  this  frequency  track 
can  continue.  However,  since  the  value  of  element  (53,3)  is  zero,  the  track 
"dies"  at  this  point. 
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mag 

This  is  the  matrix  which  contains  the  magnitudes  of  DFT.  The  number  of  rows 
of  this  matrix  is  equal  to  one  half  of  the  DFT  points,  and  the  number  of  columns 
depends  on  the  sound  signal. 

phas 

This  is  the  matrix  which  contains  the  phases  of  DFT.  It  has  the  same  size  as  the 
mag  matrix. 
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•  sig 

This  is  a  windowed  version  of  the  original  sound  signal. 

•  par 

This  vector  is  composed  of  the  essential  parameters  provided  by  the  user 
including  the  sampling  frequency,  windowed  frame  length,  overlap  width,  and 
the  number  of  DFT  points  which  were  used  in  the  synthesis  and  other  programs. 

2.  Synthesis  Function 

[sig,  synt]  =  synth(cand,  mag,  phas,  sig,  par); 

Inputs        All  input  parameters  are  generated  by  the  analysis  function. 

(see  above  description) 
Outputs     sig  windowed  signal 

synt        synthesized  signal 

B.         FREQUENCY  TRACK  EDITING  FUNCTIONS 

Incand,  nmag,  nphas,  nsig,  npar,  nsynt]  =  guifreq(cand,  mag,  phas,  sig,  par,  synt); 

This  is  the  main  program  which  needs  to  be  invoked  if  the  user  wishes  to  perform 
the  frequency  track  editing.  This  function  calls  the  associated  function  guiedit.m  which 
performs  all  GUI  and  "frequency  track  editing"  operations.  The  inputs  are  the  same 
variables  used  in  the  analysis  and  synthesis  functions  mentioned  earlier  and  the  outputs  are 
modified  versions  of  these  variables.  For  example,  ncand  is  the  resulting  candidate  matrix 
after  the  original  candidate  matrix  cand  has  been  "edited."  The  functions  called  by 
guiedit.m  are  described  below: 
edit.m  eliminates  short  frequency  tracks 

zero.m  eliminates  a  specific  range  of  frequencies 

cut.m  deletes  a  small  portion  of  signal 

paste.m  cuts  a  small  portion  of  signal  and  pastes  it  at  a  specific  time  instant 

mousedit.m       implements  the  mouse  editing  capability 


C.         SOUND  TRANSFORMATION  FUNCTIONS 

[ newsy nt|  =  guimod(cand,  mag,  phas,  sig,  par,  synt); 

This  is  the  main  program  which  calls  the  associated  MATLAB  function  guiscale.m 
in  order  to  perform  the  sound  transformations,  including  time-scale  modification  and 
frequency  scaling.  The  inputs  are  from  one  of  the  previous  two  programs,  namely,  the 
sound  analysis/synthesis  or  frequency  track  editing  programs,  and  the  output  is  the  modified 
waveform.  The  function  guiscale.m  includes  all  GUI  and  calls  function  modsynt.m  which  is 
required  to  perform  sound  transformations. 

The  user  can  also  get  on-line  help  by  typing  "help  funcname"  in  the  MATLAB 
workspace  to  look  at  the  description  about  how  to  use  the  above  functions.  Alternatively, 
one  can  simply  activate  these  interactive  GUI  tools  and  press  the  "Help"  button  to  get  more 
information. 
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